TITLE: 

TYPE: 
DESCRIPTION: 



DURATION: 



UNIVERSITY OF IIUNOIS 
DIGITAL COMPUTER LABORATORY 
ILLIAC PROGRAM LIBRARY 

Library Routine M 2k - 260 

Complete Linear Matrix Equation Solver and General 
Matrix Inversion Routine Using Drum Storage. (DOI Only) 
Complete Program (DOI built - in) 

This routine solves the matrix equation A X * B given 
A and B where 

1. A Is a non-singular square matrix of order n 

2. B ts of size (n x m). 

The magnitudes of n and m are governed by 
n + m < 164 
n 2 + 2nm + 3n < 20,^20. 
In catie of inversion (B = i) n = m and the limit is 

n < 82« 

There *#& two options on input and two options on 

output : 

Inputs 1- From tape 
2. From drum 
Output: 1. Tape output 

2, Results stored on drum. 
This routine uses the same method of solution as M 13. 
Characteristic tfcnes for calculation (input and output 
not included) are: 

Inversion: n = 20 
n = 40 
Linear equation: n = 10, m = 1 - 
n = 10, m = k - 



~ 90 sec. 
— 600 sec. 
~ 8 sec. 
~ 15 sec. 



For higher orders, the time increases approximately 

3 
as n . 

Note: Due to the large amount of drum recording necessary 
in this routine, it is not advisable to use the routine 
for low order matrices, which could be accommodated by 
M 13 or M lk which are much faster. 
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TAPE PREPARATION: Each matrix equation to be solved uses some tape input: 

either a data tape containing the matrix or a parameter 
tape, depending on the option. In all cases scaling is 
the same as in M 13. 

Input Option 1 ; When matrix A and/or B are to be read 
from tape, they appear in the form 

± a ll ± a i2 ± a l 5 — ±«m C 

±' b ll ± b I2 •••i b lm C 

± a 21 ± a 22 ± a 2n C 

± b 21 ± b 22 •••±*2m C 

i a nl ± a n2 -••± a nn C 

± b nl ± b n2 — ±* m C P 
where + a and + b are the elements of the matrices 
A and B.. C is a control character for calculations. 
P is a control character for output. 
If the problem is the solving of a linear matrix 
equation (as opposed to inversion) then the C after 
each rov is to be the sexadecimal character N. If the 
problem is the inversion of A, C is to be a J and the 
elements of B are not punched. B = I will automatically 
be constructed by the program when C = J. 
See OUTHJF for description of the P digit. 
With this option, drum Ideations 2560 to (2559 + m(m + l) 

n(n + l) * 
+ — j_- c ) vill be used by the program. 

Input Option 2 : The matrices A and/or B can be stored 

on the magnetic drum (as from previous programs) prior 

to reading in this routine. The matrices must be stored 

by row in same order as they would be punched on tape 

with a sum check after each row of each matrix. (This 

routine uses Y 1 for drum transfers.) In the case of 

inversion B need not be stored with A on the drus. 

The input tape has form 

C + n + m + dNP 
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where 

C: Control digit where Fs Solves linear equation (B f i) 

L: Inversion (B = l) 
n: Order of Beatrix A 
m: Columns in matrix B 
d: Drum location of first element of A 
N: : Sexadecimal character N 
P: Print control character (See OUTPUT ). 
Drum usage with this option is as follows : 

1. If matrix A is to be inverted (C = L), it should be 
stored above the drum location given by 

25 6l + £fe-+_il 

to allow room for forming B = I on the drum. 

2. The upper limit on drum storage is 127^9 • 
3* Working storage on the drum is from 2560 to 

(2559 + n(m + 1) + n ( n + X) . ). 

OUTPUT: The digit P determines the form of output. If 

1 < P < S 
output will be on tape by column where P determines 
number of digits to be punched for each element. Each 
column will be followed by a scaling factor and an N« 
(See M 13 for complete description of this mode). 
. If P = L the results of the computation will be placed 
consecutively on the drum by row, starting in 2560, 
a sum check following each row. After the last row the 
scaling constants for each column will be written as 
an extra row of the matrix o 
NOTE: 1. Upon completion of one problem, the program stops on a 

2^06K. 
A black switch start will start a new problem. 
2. V A stop with 
YF 010 
k2 — 
signifies a sum check st6p in YL 
3.. A stop with . 

FF Oil 
signifies a sunr check error in reading in the program. Another' 
attempt should be made to read it in. 



LOGATIOI 



10 



II 



ORDER j 

Library Routine XI- 218 
00 9K | 

Library Routine H 12 - 225 

00 48K | 

Library Routine Y 1 - 199 



00 88K 

Library Routine P 2 - 52 



00 106K 
4l 6F 

41 7F 
50 368F 
50 1L 
26 9F 
LO 250L 
36 20L 
40 8F 
L5 256L 
40 52L 
22 57L 
K5 $LL 

42 12L 
L5 3foF 
46 44l 
46 ^L 
46 i|lL 

46 a$7L 
46 $07L 
10 g()F 
42 33L 
42 X51L 
42 X75L 
42 257L 



HOPES PAGE 1 M 24 



Decimal Order Input 



Infraput 



Transfer Block of Words from the 
Memory to the Drum or from the Drum 
to the Memory 

Print (A) with or without Sign to 
n Places as Determined by a Program 
Parameter 



Clear Counters 



Call in first row 



Use Drum? 



Store I and S 1 to 



'Set y addresses 





imsm 


dra» 




■ ' MCttJttS" 
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12 


:M 26il 




1 ' ■ '■ ■■ ' ', ' " ' 






^ 


22 F 
K5 567F 


6 










42 19L 




■ 




ik 


15 30F 

46 76L 










15 


46 lift, 
46 258L 




1 






16 


10 20F 
42 80L 




>Set t addresses 

[ 






17 


42 83L 
42 86L 










18 


42 91L 
42 97L 










19 


42 259L 




' 








22 F 


13' 


J 






20 


LO 250L 
40 8F 


(5) 


\ 






21 


52 368F 
50 21L 










22 


26 9F 
L5 568F 




1 






23 


42 5F 
00 20F 




nj 

I Set Drum parameters 






24 


46 33L 
L5 369F 




( and S 1 tod 






25 


42 4F 
00 20F 




ml 






26 


46 46L 
L5 370F 




d J 






27 


L4 256L 
40 32L 










28 


L5 256L 
LO 32L 


j 


S 1 = 0? 






29 


36 36L 


J 






1 


L5 13L 


/ 

1 







M24 



LOCATIOK 


ossnt 




.fc083B r . 


PAGEE 3 


50 


14 35L 
46 50F 




/Set 7 In (112) 




31 


50 568F 
50 51L 




] 




52 
33 


26 48f 
00 F 
00 F 
F5 52L 


)4« 47* 
[27 1 

24 


1 Read row r of A 
I frcan Drum. 




3* 


1*4 5F 
40 52L 








55 


40 45L 
22 40L 




^ 




56 


50 568F 
50 56L 




1 Read row r of A 
| from tape. 




57 


26 9F 
L5 50F 




J 
\ 




58 


LO 56L 




) 




. 


10 20F 








59 


42 5F 
42 4f 




f Set n, a = n 




Mo 


22 40L 


• 


) 




41 


50 40L 
22 5L 




/ Set j addresses 




42 


F5 8F 
52 55L 
L5 256L 




} B = I? 




43 


LO 52L 
56 49L 




J S 1 ■« Ot 




44 


50 F 
50 44L 


7 


) 




45 


26 46F 
00 F 


f 55 






46 


00 F 
F5 45L 


26 


\Read roif f of B 
from Drtttn. 




47 


L4 4f 










40 52L 









M 24 



rasHar 



*9 
50 
51 

52 

53 

5* 

55 

56 

57 

58 

59 

60 

61 

62 

63 
64 

65 



1*0 *53- 
26 59L 
5oy 

50 *9L 
26 9P 
1$ 30F 
m ^9L 
10 20F 
*2 k7 
50 5L 
$6 13L 
klJ 
J9 53L 
te 53L 

10 2571* 
32 53L 
I0 257L 
IA6f 
*2 5^L 
I£25*L 
*0 7 

ia *r 

50 250L 

00 2oy 
^ 30F 

22 52L 
22 62L 
10 3F 

26 6*X 
00 207 

kS 78L 
*6 119L 
15 25% 



HOIKS 



PACK 4 



X2' 



10 5^ f 



57* 



Besd row r of B 
from tape. 



let 



Set t addresses 
(to 62L) 



Clear y 
to y + n. • 1 



Set 
B - I 



Augsent 
Avith I 
l/lO 



Set t in I 12 



a + m 



Preset Drum 





LOCATIOir 


ORDER 


I0335S PJU^ 5 




66 


40 T7L 

4o 118L 




«o - 2560 1 




61 


L5 I69L 
L4 7F 




v 




68 


lt2 84L 
42 90L 








69 


42 94L 
^2 111L 








70 


42 115L 
42 139L 




>Set x addresses 




71 


00 20F 
46 8lL 








72 


46 83L 
46 87L 








73 


46 96L 








74 


46 138L 
L3 6f 




) 




75 


l£ 7F 
56 130L 




>i>r? 




76 


50 F 
50 76L 


14* 


(call in rov i 




77 


26 48F 
00 F 


) 66 






78 


00 F 
26 79L 


64* 123* 


) 

> Reset interchange 




79 


L5 251L 


; 






46 9IL 








80 


46 92L 
L3 F 


16- 






81 


l6 f 
32 84l 


71* 


■|x| > |t|7 




82 


47 91L 


r 

L, -> 
3 






50 250L 






83 


L5 F 


72 


rCoanp ■ 










66 F 


17 


1 



K2k 



LOCATION 


j ORDER [ 
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Qk 


26 88L 








50 F 


68 


1 


85 


S3 696F 
32 87L 




V x = 0? 


86 


50 253L 




1 . 




75 F 


17' 


comp 




X 




87 


66 F 


72' 






K7 92L 




\^° 


88 


kl 5F 
Si: 532F 




Clear M box 


89 


kO 2F 




Store - k 




kl F 




Clear Sum "box 


90 


22 90L 




Waste ~\ 




L5 F 


68 r 101 


\ 


91 


40 F 


79' 82 








L5 F 


18 ,,90' 


* Interchange 




92 


kO F 


80 87 ' 






93 


50 2F 
7 J IF 




J 

, t ■ - kx 
> x ~ x - kt w 


Eliminate between 




IA F 




rows i and r 


9k 


kO F 


^/100' 


) 




95 


L3 F 
L6 5F 
36 97L 


69) 


k Store |x' I I 
1 'mx I 


96 


L7 F 


73' 101' 






kO 5F 




) \ 


97 


L5 IF 




/ Complete 






kO F 


18' 98' 


j change 




98 


F5 97L 
k2 97L 


« 


J 




99 


i*-2 91L 




^ Advance 1 




L5 9^L 




addresses j 


100 


Ik 251L 








4o 9^L 








101 


k2 90L 
46 96L 




> 






LOCATION 



ORDER 
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102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 

113 
114 
115 
116 

117 
118 

119 



LO 258L 

32 90L 

L3 5F 
36 248L 
L5 252L 
40 112L 
LL 5F 
32 108L 
L5 89L 
46 112L 
26 111L 
F5 112L 
42 112L 
L5 5? 
00 IF 
40 5F 
LL 5F 
32 107L 
50 250L 
L5 F 
10 F 
00 F 
50 250L 
40 F 
F5 H3L 
42 113L 
42 111L 
LO 259L 
32 111L 
26 117L 
JO F 
50 117T. 
26 48F 
00 F 
00 F 
L5 119L 



69' J 115 

106 '? 104' 
108J 

70 114 ' 



15 135' 






66 « 

121 



65 233 



► More elements 



) 



1 



-Singular? 



a = 
.-. ' r . 'mx 




Scale 
row 
i 



Scale to 



a > 77 

1 r. 'mx — 4- 

1 



J 



Put row i back 
on drum 



LOCATIOI 


ORDER 


ions 
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120 


10 20F 
F4 118L 








121 


4o 118L 
ko 77L 




- .-- ■■■ - 




122 


L5 119L 
LO 251L 








123 


46 119L 




) 




124 


L5 49L 
46 30F 




1 




125 


22 125L 
50 125L 




^Set y addresses 




126 


22 5L 
L5 258L 




) 




127 


46 307 
50 127L 




r""™" 




128 


26 13L 
F5 7F 




\ 




129 


42 J¥ 
26 67L 




>i -4i + 1 




130 


L3 7F 
L6 6f 




> i * r? 




131 


36 135L 




1 




132 


F5 6f 
42 6f 




/ r ->r + 1 




133 


LO 3F 
32 145L 
4l 7F 




> r = n? 

S Clear 1 and read 




134 


26 28L 

00 F 




f in next row 




135 


L5 138L 
46 117L 




\ 




1 136 


L3 6F 
32 137L 








137 


26 ll.TL 
41 5F 








* 











M24 



LOCATION 



lp-5 



1.^4 



'^>'J 



OEDKk 



158 


L3 F 




l£ 5F 


159 


32 l*OL 




L7 F 


11*0 


*0 5? 




F5 139L 


lln 


*2 139L 




00 207 


1^2 


*6 138L 




F5 7F 


1^3 


H2 ft 




ID 3F 


ltt 


IX) *F 




36 103L 


1*5 


26 138L 




81 *F 


lW 


00 20F 




*6 210L 


1*7 


L* 25IL 




00 15F 


1*8 


32 1*9L 




49 8F 


1*9 


26 150L 




in 8F 


150 


^1 5F 




L5 25*L 


151 


1*0 F 




L3 F 


152 


36 2*8L 



I£ il8L 
1*0 lb2L 
L5 119L 
1*6 l$3L 
1*1 6f 
in tf 
L5 169L 



ICOTS 
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M 2* 



7* 1*2 



70' 1*1 



8 

10' 



Bit row r 
on Brim 



1 -*i + 1 



i « at 



J 



Set punch 
and 32 



Clear k 

Fetch S Q « l/lO 
Store S 



Sing? (S < 2" 3> ) 



Set initial Drum 
addresses 

Clear counters 



1 LOCATION 


ORDER 
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156 


42 209L 




| 




L5 162L 




157 


FO 4F 
FO.. 6F 






158 


1*0 162L 




> Adjust drum 




40 193L 




addresses 


159 


L5 163L 
L4 251L 






160 


46 3-631, 


s 


i 




46 194L 


a 


; 


l6l 


50 532F 
50 l6lL 


" 


\ 


162 


26 48F 
00 F 


":i55 
158 


r""" 


±03 


Ou jc 


•154 160 


1 




> X L5 88L 




^ 


l64 


f4 6f 








42 170L 




>. Reset addresses 


165 


L4 5F 

42 167L 




1 


166 


L5 16 7L 
46 17OL 






167 


50 F 


8< 


—1 




71 F 


165' 


1 


168 


40 F 








22 17^L 




- 


l69 


00- 32JF 
§5 368F 


~v 


> Compute 5* i ir r 

■ 


170 


50 F 
74 F 


l66 ' ]l75 

l64'j 




171 


L4 f 




/ 


172 


40 F 
LL F 




}ra<§< 




32 17^L 




173 


50 25^L 




|„ _, «-. 


1 


7 J F 


11 



LOCATION 



ORDER 
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174 


26 151L 




L5 17OL 


• 175 


LO 251L 




4o 170L 


176 


42 180L 




42' 183L 


177 


46.185L 




F5 .?F 


178 


40 JF 




L5 6f 


179 


LO 7F 




32 I69L 


180 


4i 7F 




L3 F 


181 


36 248L 




1*6 f 


. iae 


36 173L 




26 183L 


183 


L5 F 




66 F; 


184 


22 184L 




SI 552F 


185 


40 F 


. 


L3 8F 


186 


32 19I1L 




L3 5F 


187 


32 194L 




L> 185L 


188 


L4 169L 




46 191.L 


169 


L5 184L 




I/i ?F 


190 


L4 6F 




42 191L 


191 


Ifr F 


_i 


40 F 

— , .._ 



176 



176 



177 



188 ' 

190' 



Advance 
addresses 
(n - r +1) times 
for row r 



Reset i 



Sing? Sii = 



Division bad? 



r a^ 



IT 



S2 = 0? 



k = 0? 



Restore row (o - r) 

of augmented matrix 

to drum *ith x , , 
r,k~l 



LQCATIOa 
192 

193 

194 

195 

196 

197 
198 

1j9 
200 
201 
202 

203 
204 
20p 
206 
20Y 
208 
209 



ORDER 

JO 532F 
50 192L 
26>8F 
00 F 
00 F 
F5 6F 
42 6F 
LO 3F 
56 197L 
22 I56L 
L3 8F 
36 209L 
L5 260L 

40 200L 

41 F 

26 200L 
L5 F 
40 F 
L5 200L 
L4 251L 
40 200L 
F5 F 

42 F 
LO 4F 
36 205L 
26 200L 
L5 206L 
L4 5F 
42 207L 
50 860F 
lO F 

40 F 

26 2I4L ' 

00 567F 

92 131F 
L5 F 



NOTES 
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158* 
160* 



) 198' 
T202 



9 
206 



'Count n rows 



S2 * 0? 
Ho 



Store co|.UBn 
k in Wif (700 to TOQfn) 
and save scaling 
factor, S. . 



156 212 



\ 



Yes 



LocAjrccar 


ORDER 
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" .' 




210 


50 F 
50 210L 


1*6 « 




211 


26 88F 
F5 209L 




/ Punch coIumb k 


212 


k-2 209L 








FO 26lL 




1 


215 


56 209L 
92 770F 




v 


21* 


F5 5F 
k2 5F 






215 


LO k¥ 




\k-*fc+l 




52 216L 




1 k fi a? 


216 


22 150L 




J 




L5 8F 




> 32 = a 


217 


52 2kQL 
L5 255L 






218 


ho 228L 

kO 2Y;l 






219 


L5 k¥ 
00 20F 




/ 


220 


k6 256L 
k6 2kjl, 






221 


L5 *F 








IA 5F 




\Set initial, addresses 


222 


00 20F 
k6 2231: 




/ for VrvB$ "Iransfer 


225 


10 20F 
IA 208L 






224 


k& 255I« 
I£ 207L 






225 


k6 25^L . 
L5 85L 






226 


46 23: L 








kl 2F 




) 


227 


50 5681- 




) 




50 22 7L 




J : m 1 



LOCATIOH 
228 

229 

2^0 

231 

252 

233 

234 

235 
236 

237 

238 

239 

240 



242 



2k 



2kk 



\ QEBER 

26 48F 
00 F 
00 P 
L5 229L 
10 20F 
F4 228L 
40 228L 
L5 229L 
LO 251L 
46 229L 
L5 F 
40 F 
JO F 
50 234L 
26 48F 
00 F 
00 F 
L5 234L 
ID 251L 
46 234L 
L5 235L 
F4 4F 
40 235L 
40 246L 

L5 233L 
L4 2^1L 
46 233L 
DO 25OL 
42.233L 
F5 2F 
42 2F 

lo 4f 

36 245L 

26 C27L 



HOTES 



PACE 14 



M 24 



} 



218 
231 

222* 232 



If 



226 24i 

224 242 

225 237' 

218 « 
239 

220 



Isolate X on 
Drua by 
row 



LOCATION 


ORDER 
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2^5 


JO 86OF 


^ 




50 2^5L 






2^6 


26 ^8f 


p9' 


Store scale 
/ constants on 




00 F 


2^7 


00 F 
22 2k&L 


220* 


drum 


2^8 


92 898F 


' F 




92 139F 




LF and CR (3) 


2^9 


2k L 
00 F 


- 


REPEAT 


250 


00 F 
00 2F 






251 


00 IF 
00 IF 






252 


10 IF 
00 IF 






253 


7L lf095F 
LL ^095F 






25*1- 


00 F 








00 1000 0000 OOOOJ \ 


255 


26 48F 
00 2560F v 




y Constants 


256 


26 48F 
00 F 






257 


K6 13L 








kl F 


11' 




258 


NO F 
L3 F 


iy 




259 


JO 250L 




1 




kO F 


19 


260 


L5 368F 
ko 696F 






26l 


12 131F 


^ " 


1 




L5 F 


12 






LOCATION 


ORDER 
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262 


1+1 F 
L5 3F 




Sum check routine 




265 


L6 f 
ho F 








263 


F5 368F 
^2 368F 








26k 


LO 375F 
32 368F 








265 


L3 F 
36 2^9L 








266 


82 ^OF 
FF 17F 








267 


22 373F 
OO F 








868 


Nl F 
L5 376F 








269 


N3 17 i +6F 
50 3893F 

26 368N 




Z check constant 




■ 




1 







RTi H/3/60 



DATE April 1, 1959 



PROGRAMMED B Y S. J. Penny 
APPROVED BY Jo H. Snyder 



nj 



